Shot noise in the chaotic-to-regular crossover regime 
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We investigate the shot noise for phase-coherent quantum transport in the chaotic-to-regular 
crossover regime. Employing the Modular Recursive Green's Function Method for both ballistic 
and disordered two-dimensional cavities we find the Fano factor and the transmission eigenvalue 
distribution for regular systems to be surprisingly similar to those for chaotic systems. We argue 
that in the case of regular dynamics in the cavity, diffraction at the lead openings is the dominant 
source of shot noise. We also explore the onset of the crossover from quantum to classical transport 
and develop a quasi-classical transport model for shot noise suppression which agrees with the 
numerical quantum data. 
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The significance of noise induced by the discreteness 
of the electron charge ( "shot noise" ) first investigated al- 
most a century ago |l| has resurfaced in the field of meso- 
scopic physics 2, 3]. Shot noise carries information about 
the crossover from a deterministic (classical) particle pic- 
ture of electron motion to a probabilistic (quantum) de- 
scription, where electrons behave as matter waves. The 
uncertainty inherent in a quantum picture gives rise to 
noisy transport. In conductance through quantum dots 
the correlations between electrons in the Fermi sea lead 
to a suppression of shot noise S relative to the Poissonian 
value of uncorrelated electrons Sp 4] . The reduction rel- 
ative to the completely random limit is customarily ex- 
pressed in terms of the Fano factor F — S/Sp. 
Most investigations to date have focused on quantum 
dots whose classical dynamics is fully chaotic 0, H, H, 
H, EJ El [n|. In this limit, random matrix theory 
(RMT) yj predicts a universal value for the Fano factor, 
F = 1/4. The applicability of this RMT result requires, 
in addition to the underlying chaotic dynamics, dwell 
times in the open cavity tjj which are sufficiently long 
compared to the Ehrenfest time te- The latter estimates 
the time for the initially localized quantum wavepackets 
to spread all over the width d of the cavity (typically 
d « \f~A with A area of the dot) due to the divergence of 
classical chaotic trajectories. It can be estimated as [l2^ 



T E = A- 1 Hd/\ F ), 



(1) 



where A is the Lyapunov exponent (A > for a chaotic 
cavity), and Xp is the de Broglie wavelength associated 
with the wavenumber at the Fermi surface kp. The 
limit te/td <C 1 corresponds to the quantum (or RMT) 
regime and te I to 1 corresponds to the classical limit 
for which F = is expected. For ballistic cavities in the 
crossover between these two regimes a simple conjecture 
for F was put forward 0, 



F=l/4 exp(-r E /T D ) 



(2) 



alternative crossover behavior, 



F = l/4 {1+tq/t d )- 



(3) 



For cavities with a short-ranged disorder potential, an 



was proposed 0, 0] , where tq is a characteristic scat- 
tering time within which the wavepacket is scattered into 
random direction and thus explores the entire cavity, tq 
and te are closely related to another as both denote the 
characteristic time scale for spreading of the wavepacket 
by chaotic scattering either at the boundary (te) or the 
interior (tq) of the cavity. Moreover, for short ranged dis- 
order with a correlation length lc < Xp , tq incorporates, 
just as Tp, quantum effects and depends on an effective 
H e s of the system. The crossover from the chaotic to the 
regular regime is therefore predicted to be controlled by 
a single ratio te /to or tq /to which will be a function 
of both the size of quantum effects (fi e fi) and the mean 
rate of irregular (chaotic) scattering, (A). The chaotic- 
to-regular crossover corresponds to the limit (A) — > 0, 
while the quantum-to-classical limit involves ft G ff — > 0. 
An obvious test of different predictions would be to sim- 
ulate phase-coherent scattering processes on a computer. 
However, conventional numerical techniques suffer from 
a slow convergence rate for large te and, in particular, 
for small fi e ff- To circumvent this difficulty an "open" 
dynamical kicked rotator model was recently success- 
fully used to mimic chaotic scattering in a ID-system 
[(1E3- Experimental tests have, so far, been limited to 
the regime of tq/td < 1 ^3 ■ 

One open question not yet well understood is the be- 
havior of shot noise for motion in a regular rather than 
chaotic cavity, i.e. in the limit A — > 0. For a mixed system 
lower values of F have been observed |r| suggesting that 
for regular systems F may vanish. Taken at face value, 
Eq. J2J) yields F = for the case of te ~ * 00 or A — > at 
fixed value of h c g. This result would correspond to a com- 
plete suppression of shot noise. On the other hand, the 
notion of noisy transport as a result of the wave nature 
of electrons would suggest that noise should be inherent 
in phase-coherent transport through mesoscopic struc- 
tures irrespective of the chaoticity of the underlying clas- 
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FIG. 1: Rectangular quantum billiard with tunable shut- 
ters and tunable disorder potential (gray shaded area). Elec- 
trons are injected from the left into the cavity region of size 
A = 2d 2 , width 2d and height d. Tuning the opening of the 
shutters w and the strength of the disorder potential Vo, the 
onset of the crossover from quantum-to-classical and chaotic- 
to-regular scattering can be investigated, respectively. 



sical dynamics. In the present letter we analyze a model 
for transport through a cavity that allows to investigate 
the crossover regime from chaotic to regular dynamics, 
i.e. A — > 0, and the onset of the crossover from quantum 
to classical dynamics. In our model, the dwell time and 
the mean rate of chaotic spreading of the wavepacket, 
(A), can be tuned independently. Surprisingly, we find 
that the Fano factor and the distribution of transmission 
eigenvalues for regular systems resemble closely that of 
chaotic systems. We interpret this observation in terms 
of the ubiquity and dominance of diffractive scattering in 
phase-coherent cavities. 

We choose first a scattering geometry (Fig. which 
consists of a rectangular cavity to which two leads of 
width d are attached via tunable shutters with an open- 
ing width w. To reduce direct transport between the 
shutter openings, they are placed on the bottom and top 
end of the leads, respectively. The cavity region of width 
d and length 2d contains a disorder potential V char- 
acterized by its mean value (V) = 0, and the correla- 
tion function (V(x)V(x + a)) = (V 2 ) exp(-a/7 c ) 03 
The correlation length lc is typically a small fraction of 
the Fermi wavelength lc/^F ~ 0.12 and the potential 
strength V = ^(V 2 ) is weak, V /E F < 0.1 . The depen- 
dence of td on the shutter opening was determined by 
a Monte-Carlo sampling of ~ 10 5 classical trajectories, 
td ~ 2.66/(wkp), i.e. approximately independent of Vq. 
In the limit of vanishing disorder potential (Vo — > 0) the 
motion inside the cavity becomes completely regular. 
Our quantum calculation proceeds within the frame- 
work of the modular recursive Green's function method 
(MRGM) T7] which allows to treat two-dimensional 
quantum dots with relatively small Xp (or small H e g). 
The MRGM requires, however, separable modules. We 
construct the latter by decomposing the cavity region 
into two square modules for each of which we choose a 
separable random potential V(x, y) — V\{x) + V^y). In 
order to destroy the unwanted separability and to ensure 
chaotic dynamics, we build up the cavity by combining 



two identical modules, however, rotated by 180° relative 
to each other (see Fig. ^ for illustration). As the Fermi 
wavenumber and the width of the leads d is independent 
of the shutter openings w, this device allows to study 
two different crossovers of shot noise: (1) By changing 
the opening of the shutters w the dwell time in the cav- 
ity td can be tuned without changing the Fermi energy 
Ep. With the increase of w the onset of the quantum- 
to-classical crossover can be probed, with the sidcffcct, 
however, that not only H e g is reduced but simultaneously 
the classical phase space structure is altered. This closely 
resembles the parameter tuning in the experiment |l4j . 
(2) By tuning the strength of the disorder potential, the 
rate of chaotic spreading, (A), and thus Tp (or tq, used 
in the following interchangeably) is varied. 
We evaluate the transmission amplitudes t mn for an elec- 
tron injected from the left by projecting the Green's 
function at the Fermi energy G(Ep) onto all modes 
m,n £ [1, . • ., N] in the in- and outgoing lead, respec- 
tively. The Fano factor F is then calculated from the 
A^-dimensional transmission matrices t 

(Trttt(l-ttt)) _ (E^TUI-TJ) 

with T n being the eigenvalues of t't. The brackets (. . .) 
indicate that we average over 150 equidistant points in 
the wavenumber-range kp € [40.1,40.85] x n/d, where 
40 transverse lead modes are open. Note that this mode 
number is higher than in previous studies. We choose a 
nearest-neighbor spacing Ax = Ay in the Cartesian dis- 
cretization grid such that the Fermi wavelength is well 
resolved by a large number of grid points, Xp w 32 Ax. 
With these settings we have a total number of ~ 8.5 x 10 5 
grid points in the interior of the cavity. 
Figure |3 displays the Fano factor as a function of the in- 
verse dwell time r^ 1 (Fig.Et) and the inverse Ehrenfest 
time t^ 1 (Fig. Eh)- For r^ 1 — > (i.e. large dwell times) 
F approaches the universal value 1 /4 irrespective of the 
strength of the disorder potential Vq, while for shorter 
dwell times F falls off gradually (Fig. . The steepness 
of this decrease is clearly dependent on V and thus on 
the mean scattering rate (A). Most striking is the fea- 
ture that for Vo — > but long dwell times the shot noise 
reaches the RMT value even though the dynamics is now 
entirely regular (see Fig. (2Ji). In Fig. 03 this feature is 
reflected by the fact that the Fano factor F does not de- 
cay to zero even as Vo — > (i.e. r^ 1 — > 0). There is 
no compelling a priori reason why the RMT prediction 
should be applicable to regular systems. This observa- 
tion suggests that the conjectures [Eq. © or Eq. 
require a modification to properly account for the shot 
noise in the regular limit. We argue that the key point is 
the wavepacket diffraction at the cavity openings which 
has to be incorporated in the theoretical description of 
shot noise. Note that this feature is inherent in quantum 
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FIG. 2: (a) Fano factor F in the quantum-to- 
classical crossover regime for the cavity depicted 
in Fig. Curves shown correspond to differ- 

ent strengths of the disorder potential: Vo/Ef = 
0.1(h), 0.07(d), 0.05 (•), 0.03 (o), 0.015(a), 0(a). 
(b) Regular-to-chaotic crossover regime of F. Curves cor- 
responding to different shutter opening ratios are shown: 
w/d = 0.55(h), 0.5(d), 0.4(»), 0.3 (o), 0.2(a), 0.1(a). In- 
sets of (a) and (b) depict the fit parameter free prediction with 
diffractive corrections [Eq. JjJ]. (c) Distribution of transmis- 
sion eigenvalues p(T) in the quantum-to-classical crossover: 
Vo/E F = 0(T),0(V) and w/d = 15(T),50(v). (d) p(T) in 
the chaotic-to-regular crossover: Vo/Ef = 0.1(4), 0(0) and 
iu/d = 50f»,50 (0). 



estimate for Tp (or tq). We perform a quasi-classical 
Monte-Carlo transport simulation [2(| in which we fol- 
low an ensemble of classical trajectories subject to a ran- 
dom Poissonian scattering process. For the latter we 
calculate the transport mean free path (rg ■ Vp) and 
the differential scattering probability (P(6) ~ da/d9) 
in first Born approximation, thus taking into account 
quantum diffractive scattering (for lc • kp < 1) along 
the lines of Refs. 0,0]. In this simulation, spreading 
results from both the injection of classical trajectories 
with an initial angular distribution given by Fraunhofer 
diffraction at the shutter opening and multiple scatter- 
ing inside the cavity. This allows to identify a modified 
Ehrcnfcst time Tp as the time it takes for the ensem- 
ble to acquire a mean spread of the order of the cavity 
width d under the influence of both disorder (or chaotic) 
scattering inside the cavity as well as diffraction at the 
shutter openings. At te the mean separation reaches 
Ar(TE) = ({r— (O) 2 ) 1 ^ 2 = d. Already the inclusion of 
Fraunhofer scattering alone leads to a drastic reduction 
of the Ehrenfest time: Injected electrons spread much 
faster over the whole cavity than according to the esti- 
mate in Eq. Q . An important consequence of this mod- 
ification is that the parameter rg is now of comparable 
magnitude or even smaller than the time To beyond which 
universal behavior for P(t) sets in [P(t) oc exp (— tjro) 
for chaotic and P(t) oc {t/To)~^ for regular dots]. In this 
regime system-specific deviations of the dwell time distri- 
bution from a universal decay law are more pronounced 
than characteristic mean differences between regular and 
chaotic cavities 0, 0, E| . For our estimate of the Fano 
factor F we therefore take into account the exact dwell 
time distribution P(t), resulting from the quasi-classical 
Monte-Carlo trajectory simulation for the particular sys- 
tem we study. Following the analysis jfj based on the 
contribution of noiseless transmission channels, we find 



F= 1/4 



1 - 



P(t) dt 



= 1/4/ P(t)dt. (5) 



transport and independent of the underlying regular or 
chaotic dynamics |TjJ . Scattering due to chaotic dynam- 
ics, which lies at the core of RMT, certainly leads to 
wavepacket spreading but does not constitute the only 
or, in general, dominant source. 

Diffraction at the cavity openings has been studied in 
detail 0,0] and can be described by a standard Fraun- 
hofer diffraction analysis for electrons which enter or 
leave the cavity ■ An estimate for wavepacket spread- 
ing can be found by considering the characteristic angu- 
lar injection patterns of the transverse modes in these 
openings. Averaging the individual patterns over the 
total number of modes in the shutter openings, M — 
int(fcgu>/7r), we obtain for the variance of the injection 
angle 9 the dependence: (^(6 2 ) - (6>) 2 ) w 0.5 x M~ - 5 . 
This result now enters our considerations on a modified 



Note that Eq. (JSJ is applicable to chaotic and regular sys- 
tems and is expected to be valid irrespective of whether 
the origin of spreading is ballistic scattering at the bound- 
ary or diffractive scattering inside of the cavity. 
It is now tempting to probe the convergence towards the 
RMT limit for te/td < 1 also for the distribution of 
transmission eigenvalues p(T) and for other regular and 
chaotic structures that feature only ballistic motion. We 
therefore investigate transport across the semi-circular, 
rectangular and semi-stadium shaped cavity which are 
prototypical for regular and chaotic ballistic motion, re- 
spectively (see insets of Fig.|3jD,c,d). Unlike in the case of 
the geometry of Fig.^ where the number of open modes 
was fixed, transport coefficients were calculated here for a 
wide range of hp- The dynamics of these cavities is clas- 
sically scaling invariant. An increase of kp thus probes 
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FIG. 3: (a) Fano factor F for scattering systems with chaotic 
(stadium) and regular (circle, rectangle) classical dynamics 
as a function of kp, in units of the number of open modes 
N — mt(kpb/n) (each data point is averaged over 200 equidis- 
tant fcF-values). The inset shows the fit parameter free 
prediction with diffractive corrections [see Eq. (b)-(d) 
Distribution of transmission eigenvalues p(T) averaged over 
the Arrange depicted in (a). Remarkably the obtained val- 
ues for p(T) correspond very closely to the RMT-prediction 
P 

insets show the cavity geometries 



' (T), irrespective of the chaoticity in the cavity. The 



directly the dependence of Tp/rp on fi, e ff, while keeping 
the classical phase space and A fixed. Unlike the Tp — > 
crossover performed with the help of the tunable shutters 
as in the experiment |l4| and in the simulation above, the 
kp — > oo limit would correspond to the "pure" quantum- 
to-classical crossover as fi. e ff — ► 0. 

Employing the same high density of grid points we reach 
up to 15 flux-carrying modes in the semi-circular and 
rectangular billiard and 7 open lead modes in the case 
of the semi-stadium . We therefore can observe only 
the onset of the quantum-to-classical crossover for these 
three cavities. Their area A = (4 + 7r)/2 and lead width 
b = 0.125 are identical. Given this small leadwidth, 
i.e. bj 1 \J~A = 0.066, our results agree very well with those 
shown in Fig. [3] for small shutter openings: at all ener- 
gies the Fano factor is close to 1/4, irrespective of the 
chaoticity in the cavity. However, for a higher number 
of open modes (i.e. higher kp) diffraction becomes less 
important since te becomes larger and a decrease in the 
Fano factor ensues, signifying the onset of the quantum- 
to-classical crossover. Our numerical data (Fig.|3Ji) agree 
indeed with this prediction. Moreover, the quantitative 
behavior is well described by Eq. (see inset of Fig. EH). 
The RMT prediction for the distribution of trans- 



l/[wy/T(l — T)]. Taking our original scattering geome- 
try (Fig.^) with a setting of wide shutter openings (short 
Tp, or "classical") and no disorder potential ("regular") 
we find a curve which clearly deviates from this RMT- 
prediction (see Fig. |2H,d). Remarkably, the RMT-limit 
is, however, restored equally well by either (1) further 
closing the shutters (longer Tp or "quantum") or (2) 
by increasing the disorder potential (increase of (A) or 
"chaotic"). In line with this observation the stadium- 
shaped, circular and rectangular cavity all fulfill the 
p RMT (T) remarkably well (see Fig. Et>,c,d) due to their 
small lead openings (i.e. long dwell times). This suggests 
that for Tp/rp, <C 1 not just F but also the distribution 
function p(T) approaches the "universal" RMT limit. 
To summarize, we have numerically verified the behav- 
ior of the Fano factor F in a realistic scattering system 
with a tunable disorder potential and tunable shutters. 
We find that diffractive quantum scattering is sufficient 
to establish the RMT eigenvalue distribution p RMT (T), 
irrespective of regular or chaotic dynamics. The chaotic- 
to-regular and quantum-to-classical crossover in F can 
be estimated by a generalization of a previously proposed 
dependence [jj on the Ehrenfest time fp [Eq. J5J)], pro- 
vided that the definition of the Ehrenfest time is properly 
modified to include diffraction. 
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